Amendments to the Claims: 

This listing of claims will rqjlace all prior versions, and listings, of claims in the 
application. 

Listing of Claims: 

Claim 1 (original): Method for estimating the seismic illumination fold | (x, p) in the 
migrated 3D domain at at least one image point x, for at least one dip of vector p. 

Wherein the illxmiination fold I (x,prs,r) for each (source s, receiver r) pair in the seismic 
survey is estimated, by applying the following steps: 

- determination of the reflection travel time tr (Xr (p); s, r) from the source "s to the 
specular reflection point x, on the plane reflector passing through the image point Ti and 
perpendicular to the dip vector p and then returning to the reflector r; 

starting from the diffraction travel time ta (x; s, r) from the source 1 to the said image 
point x and then returning to the reflector f ; 

- incrementing the said illumination fold l(x.p;s,0 related to the said (source s, receiver 
T) pair as a frmction of the difference between the diffraction travel time U (x;s^ and the 
reflection travel time tr (Xr (p); s. 0- 

Claim 2 (original). Method according to claim 1, comprising the step of summating each 
of the said illumination folds I(x,p;s,r) related to a (source s, receiver r) pair so as to determine the 

total illumination fold I(x,p) = X i(x,P;s,f). 

sp 

Claim 3 (currently amended): Method according to claim K ono of the pr e ceding oloims, 
wherein, during the incrementing step, the illumination fold I Oc,p,s,f) is incremented using an 
increment frmction I (ta, tr; s,r) according to I(x,^ = I(x,^ + i (td, tr; s,f), the said increment 
frmction taking account of the difference between the diffraction travel time ta (x; s, T) and the 
reflection travel time Lr (xr (p); s, r). 

Claim 4 (original): Method according to claim 3, wherein the increment frmction i is a 
frmction of the seismic wavelet s(t). 

Claim 5 (original): Method according to claim 4, wherein the increment frmction i is 
expressed as a frmction of the derivative of the seismic wavelet s(t) according to: 

I(td, tr;l,T)-s(td(x; s, r).tr(xr(p); s,T) 
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Claim 6 (original): Method according to claim 4, wherein the increment function i is 
expressed as a function of the derivative s (t) of the seismic wavelet s (t) with respect to time 
according to: 

i(td, tr; s, r)-s(td(x; s, r) - Lr(Xr( p ); s, r). 

Claim 7 (original): Method according to claim 3. a ny one of claims 3 to 6. in which an a 
priori correction w (x,s,T) of the illumination fold is taken into account by migration, comprising 
the step of incrementing the illumination fold I (x, p;l,T) related to a (source T, receiver T) pair 
by i (td, tr; T, T) . w (x; s, f ). 

Claim 8 (currently amended): Method according to claim L any one of th e preceding 
claims, wherein the determination step includes the second order Taylor series development of 
the diffraction travel time (x; s, r) around the image point x: 

td(x; s, T) = td(x; s,7) + (Vxtd(x; s, r) )^.(xrx) + 2(Xr-x)^ Ax, xta (x; s^ 

Claim 9 (original): Method according to claim 8, wherein the specular reflection point Xr 
(p) is determined along the length of the said reflector such that the diffraction travel time at the 
said specular reflection point Xr (p) is stationary, according to the equation: 

Fa (7 X td (x; s, r) + Ax, xtd (x; s, r) . ( Xr ( p ) - x )) = o. 

Claim 10 (currently amended): Method according to claim 8. any on e of claims 8 or 9. 
wherein the specular reflection point Xr and the reflection travel time t, (Xr ( p ); s, r) are 
determined according to the following expressions: 

Xr(p) - x-M.F-'.b 

tr(xr(p);s,r) = td(x;s,r) -)2 -b'. F"'.b 

where: 

- M is a (3 X 2) matrix described by two vectors extending along the length of the 
reflector, and therefore perpendicular to the dip vector p; 

- b is a (2 x 1) vector of first order derivatives of the diffraction travel time along the 
reflection plane: b = W . (7 xtd); 

- F is a (2 X 2) matrix of second order derivatives of the diffraction travel time along the 
reflection plane: F = . (A x, xtd) • M. 

Claim 1 1 (original): Method according to claim 10, any one of the pr e ceding claims. 
wherein the determination step uses isochronic migration maps td (x; s, r) specified for each 
(source s, receiver?) pair involved in the migration at each image point X in the migrated 3D 
domain. 
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Claim 12 (original): Method according to claim IL a ny one of the preceding claims, 
wherein the seismic illumination fold I (x, p) in the migrated 3D domain is estimated during the 
Kirchoff summation migration of seismic data recorded during the 3D seismic prospecting. 

Claim 13 (currently amended): Method for correction of seismic data amplitudes 
recorded during 3D seismic prospecting in order to compensate for the effect of non-uniform 
illumination of sub soil reflectors, comprising the steps of: 

- estimating the illumination fold I (x, p) using the method according to claim 1. any one 
of claims 1 to 12 , 

- using the inverse I (x, p) of the said ratio as a weighting factor to be applied to each 
of the said seismic data amplitudes. 

Claim 14 (currently amended): Method for selection of an acquisition geometry among a 
plurality of acquisition geometries as a function of the target of 3D seismic prospecting, 
comprising the steps of: 

- determining the illumination fold I (5r, p) by the method according to claim 1. any one 
of claims 1 to 12, for each of the acquisition geometries considered, 

- selecting the acquisition geometry providing the optimum illumination fold as a 
function of the target. 

Claim 15 (new): Method according to claim 2, wherein, during the incrementing step, 
the illumination fold I (x,p,s,r) is incremented using an increment function i (ta, tr; s,r) according 
to I(x,p) = I(x.p) + i (td, tr; s,r), the said increment function taking account of the difference 
between the diffraction travel time td (x; s, r) and the reflection travel time tr (Xr (p); s, r). 

Claim 16 (new): Method according to claim 15, in which an a priori correction w (x,s,r) 
of the illumination fold is taken into account by migration, comprising the step of incrementing 
the illumination fold L (x, p; s, T) related to a (source s, receiver r) pair by i (td, t,; s, T) . w (x; "s, 
^. 

Claim 17 (new): Method according to claim 16, wherein the determination step includes 
the second order Taylor series development of the diffraction travel time (x; s, r) around the 
image point x: 

td(x; s, T) = td(x; s, r) + (Vxtd(x;s,r))'^.(xrjO + 2(Xr-x)^ Ax,xtd(x^ 
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Claim 18 (new): Method according to claim 17, wherein the specular reflection point Xr 
(p) is determined along the length of the said reflector such that the diffraction travel time at the 
said specular reflection point Xr (p) is stationary, according to the equation: 

pW (Vxtd (x; s, r) + Ax, xta (x; s, f) . ( Xr ( p ) - x )) = a 

Claim 19 (new): Method according to claim 18, wherein the specular reflection point Xr 
and the reflection travel time tr (xr ( p ); s, T) are determined according to the following 
expressions: 

Xr(p) - x-M.F-'.b 

tr0^r(p);s,i) = td(x;s,r) -)2 b'. F"^ .b 
where: 

- M is a (3 X 2) matrix described by two vectors extending along the length of the 
reflector, and therefore perpendicular to the dip vector p; 

-b is a (2 X 1) vector of first order derivatives of the diffraction travel time along the 
reflection plane: b = M^ . (V xtd); 

• F is a (2 X 2) matrix of second order derivatives of the diffraction travel time along the 
reflection plane: F = M^ . (/\ x, xtd) • M. 

Claim 20 (new): Method according to claim 19, wherein the determination step uses 
isochronic migration maps ta (x;"s, T) specified for each (source ?, receiver T) pair involved in the 
migration at each image point x in the migrated 3D domain. 

Claim 21 (new): Method according to claim 20, wherein the seismic illumination fold I 
(x, p) in the migrated 3D domain is estimated during the Kirchoff sunmiation migration of 
seismic data recorded during the 3D seismic prospecting. 
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